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Abstract 



By applying the supersymmetric approach we rigorously prove smooth- 
ness of the averaged density of states for a three dimensional random band 
matrix ensemble, in the limit of infinite volume and fixed band width. We 
also prove that the resulting expression for the density of states coincides 
with the Wigner semicircle with a precision 1 / , for W large but finite. 

1 Introduction 

Random Matrix Theory (RMT) has proved to be relevant in the study of several 
physical models. It was initially applied to the study of resonance spectra of 
complex nuclei and later to the study of the quantum properties of weakly 
disordered conductors, and the spectral properties of quantum systems which 
are chaotic in their classical limit |jl||2|. RMT also appears in other fields, such 
as statistics, number theory and random permutations. See for example [DQl^ 
for recent developments. 

In this article we study the density of states for a class of Hermitian random 
matrices ifj, whose elements are Gaussian with mean zero and covariance 



In the classical case of GUE, Gaussian unitary ensembles, the indices i and j 
range from 1 to iV and = 1/iV. For this case the density of states (DOS) is 
given by Wigner's famous semicircle law 



in the limit N '\ oo. Our analysis will focus on band random matrices for which 
the indices i, j range over a box A H and Jij is small when |i — j | is larger 
then some fixed band width W. As we let A t the spectral properties of 
such matrices should be quite similar to that of a random Schrodinger operator 
on a lattice given by 



{HijHki) — Sjk Sii Ji 



(1.1) 
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where V{j) are independent random variables and A is the discrete Laplacian. 
For example in one dimension the spectra of the random Schrodinger and the 
random band matrix are pure point with exponentially decaying eigenfunctions 

|7j. This is called localization. In two or more dimensions, localization also 
holds for energies outside some interval depending on d, A, W. Thus band 
random matrices are a way of interpolating between classical random matrix 
ensembles (GUE or GOE) and random Schrodinger. 

The goal of this article is to obtain detailed information about the density 
of states for a special class of random band matrices in 3 dimensions. We shall 
consider energies at which extended rather than localized states are expected. 
More precisely let z, j G A n Z^, A a set of cubes of side W, and define 



J^i ■■= -^TTTITTT - 4^W^7TTrae - (1.4) 



where A is the Laplacian with periodic boundary conditions in the volume A and 
W is large but fixed. Our estimates are valid uniformly in the size of A C Z^. 
The average density of states is given by 

(pA(i?))=-ilimIm (( ^^^ „ ) ) ■ (1.5) 

Note that as A "f Z^, the density of states p{E) does not depend on the config- 
uration with probability one. The derivative of pa{E) is 



d 

dE 



(pa(^)) = -limlm R(E + ie;0,x) (1.6) 



where 



Note that for x 7^ 

because of the symmetry Hij —Hij. 

Our main result is that for large W and E inside the interval [—2, 2], {p\{E)) 



equals the Wigner semicircle distribution (1.2) plus corrections of order . 
Moreover R[x) decays exponentially fast and piE) is smooth in E. These results 
hold for fixed W and are uniform in e as e | and in the volume as A t Z,^ . 
See Theorem 2.1 for a precise statement. When \E\ > 2 + 0{W~^) we expect 
that p{E) is smaller than any power of and that localization holds. 



For random Schrodinger operators given by (1.3) we have the classic bound 
by Wegner, p{E) < const A^^ for small A. This estimate is far from optimal since 
for small A we expect the density of states to approach that of the Laplacian. 
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Unfortunately there are no uniform bounds on p{E) as A —> or estimates on 
the smoothness of p{E) unless either the distribution of V is Cauchy (in which 
case the density of states can be explicitly computed) or E lies in an interval 
for which localized states are proved to exist. 

Note that for both the random Schrodinger and the random band matrix 
ensembles it is conjectured that for d = 3 



1 



E + ie- H 



Dp- 



^(^^ : (1.9) 



where D is the diffusion constant. Here E must be inside [0, 4c?] for the case of 
random Schrodinger or inside [—2, 2] for our band matrix ensemble, and both 
and A are small. This paper does not address this important conjecture. 
Instead we are using the phase oscillations of the Green's functions to obtain 
exponential decay for R(x). 

To establish our results on Green's functions we use the supersymmetric 
formalism of K. Efetov [|ll| which has its roots in earlier work by Weg- 
ner |^ . We recommend the survey article of Mirlin ||l| and also the paper 
of Mirlin and Fyodorov |jl^ which studies random band matrices in 1 dimen- 
sion. In the mathematics literature, A. Klein studied the density of states using 
supersymmetric methods |p3| but only at energies where localization holds. 

The supersymmetric method enables one to explicitly average the Green's 
function over the randomness. This technique involves the use of both real and 
anticommuting variables. However when we perform our estimates all anticom- 
muting variables are integrated out so that the resulting integrals is just over 
real variables. As a result of this averaging, the problem is converted into a 
problem in statistical mechanics whose action has approximately the form 

where the potential [/ is a function of the field (j>{j) and has two saddle points. 
In some respects this problem looks like a double well cj)'^ interaction. A more 
careful analysis of the integral over (p shows that one saddle dominates and it 
yields the Wigner semicircle distribution. The second saddle is suppressed by 
a determinant as we shall explain later. The large parameter W ensures that 
the integral is governed by the saddle and its Gaussian fluctuations. There are 
similar integrals which appear for random Schrodinger operators, however the 
path integral is much more oscillatory and we can not yet control them unless 
there are long range correlations in the V{j). 

The average oi \ {E + ie — E[)q^\^ can also be calculated with the supersym- 
metric formalism but the statistical mechanics is now more complicated. Instead 
of two saddle points there is a non compact saddle manifold and fluctuations 



have massless modes which are responsible for the power law (1.9). 

The remainder of this article is organized as follows. In Sec. 2 we give a 
precise statement of our results. In Sec. 3 we use the supersymmetric formalism 
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to convert averages of the Green's function to a model in statistical mechanics. 
The advantage of this representation is that for large W we see that the integral 
is dominated by two saddle points. These saddle points and their Hessians are 
discussed in Sec. 4. The following section is devoted to obtaining our results in 
a box A C of side W. In the last section we show that the analysis in the 
box can be extended to using a variant of the cluster expansion. 

Notation As in the paper we will need to insert many constants in the dif- 
ferent bounds we will denote by K any large positive constant, independent 
from W and A, and by c any small positive constant independent from W and 
A. These constants need not be the same in different estimates. Also we will 
sometimes use the symbol < to indicate that there is a constant factor K on 
the right side of the inequality (< stands for < K) without writing K explicitly. 

Acknovkrledgments. We thank Rowan Killip for many discussions and sug- 
gestions related to this paper. These discussions lead us to improve the proof, 
in particular by introducing Brascamp-Lieb inequalities. 



2 Model 

As we said in the introduction, we consider the set Ti. of Hermitian matrices 



H with entries i,jGAcZ'^,d>0. From (LI) we see that the probability 
density is 

^(^) - n n -/^^ (2-1) 



where < is an order relation on A and J is defined in (1.4). With these definitions 
is a set of hermitian random band matrices with band width W . Note that 
in d = 1 for A = [1, TV] we have |A| = N and = N~'^ exp[-|i- j|/7V] which is 
very close to GUE. For any function of H F{H) we define the average {F{H)) 
as 

{F{H)) ^ j dH P{H)F{H) (2.2) 
We study the averaged density of states pa{E): 

where "Im" indicates the imaginary part and is the retarded Green's func- 
tion: 

G+ := ^- . (2.4) 

E + ie-H ^ ' 
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In the following we restrict to d = 3 and we consider energies well inside the 
spectrum. For technical reasons we also avoid the energy E — 0. Therefore we 
consider only energies in the interval 

I =: {E : \E\ < 1.8 and \E\ > T]} (2.5) 

with 7/ > 0. We assume that our region A is a union of cubes in of side W. 
The paper is devoted to the proof of the following theorem 

Theorem 1 For d = 3, there exists a value Wq such that for all W >Wo the 
averaged density of states pa{E) is smooth in E, in the interval I, uniformly in 
W and A (hence also in the limit A | ): 

\dlPA{E)\ < C„ yn<no{W) (2.6) 

where limv^Too no{W) — oo. Moreover pk{E) is the semicircle law with a preci- 
sion 



\pa{E) - Psc{E)\ < (2.7) 



where psc{E) is the semicircle law psc defined in (l.i). Note that the first 
equation ( |g.6[ j means in particular that for x ^ Q 

i?(x) = |(G+ G+)^^J < -^e-^^ (2.8) 



Outline of the paper In Sec. 5 we establish Theorem 1 on a cube A of 
side W . We use the supersymmetric formalism to write {G^q) as a functional 
integral where a saddle point analysis can be performed. Actually there are 
two saddles. For d = 3 one saddle is suppressed by a factor e~^ (note that 
this is not true for d < 2). The fluctuations around the saddle are controlled 
using small probability arguments while the integral near the dominant saddle 
is estimated by a Brascamp-Lieb inequality . 

Sect. 6 is devoted to the cluster expansion which enables us to analyze the 
limit A ] . The cluster expansion expresses (GJq) as a sum over finite volume 
contributions Y C which are again unions of cubes of side W . We show that 
large Y terms give small contributions to (Ggo)- This expansion also enables us 
to prove the bound on R{x). 



3 Supersymmetric approach 

In this section we shall use the algebraic formalism of supersymmetry to express 
our average Green's function in terms of a functional integral, which, apart form 
a determinant, is local. Let J be given by (|l.4|). 
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Lemma 1 The averaged Green's function can be written as 



(Gqq) = I dadb cxp 
1 

{Ee - ao) 



n 

iGA 



[Ee - lb,) 

{E, - a,) 



dei[J-^ ~ F{a,b) ^ F'{aM\ 



(3.1) 



where a and b are vectors whose components ai, b^ (i = 1,...,|A|J are real 
variables and and b^ are the corresponding transposed vectors. We defined 
= E + ie and F{a, b) and F'{aQ, fop) o.f'^ matrices with elements 



F{a, b). 



1 



[E, - a,){E, - ib,) 



(3.2) 



{E, ~ ao){E, - tbo) 



(3.3) 



Note that each Oi has a pole at Oi — E + is while bi has no singularity ( as it 
appears only in the numerator) . This expression is then well defined only for 
£ > 0. 

By the same technique we obtain a similar formulas for {G^j G^q) and in 
general for {G^^^ G+j^ •■• o)- 

Remarks Note that if we omit the observable, that is we omit {E^ — uq)^^ 
and F' in (O), we are actuaUy computing (1) = 1, thus 



1 



da db exp 



n 

.ieA 



{Ee - ih) 
[Ee - a.) 



det[J-i - F(a,6)] . 

(3.4) 



Proof Note that the Green's function can be written as a functional integral: 
G+ = y" dS*dS exp [iS+{E, - H)S] Act[-i{E, - H)] SkS^ (3.5) 



where the determinant is the normalization factor and we defined 
/ ^1 



S ■ 



{si.- 



|A| 



(3.6) 



A| 



and 5*1,... 5*1 A| are complex bosonic fields. In order to insert all the H dependence 
in the argument of the exponential we introduce integrals over fermionic fields: 



Aet[-i{Ee -H)]= J dx*dx exp [tx^{E, - H)x\ 



(3.7) 
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where 



V X|A| 



>X|A| 



(3.8) 



and Xij---X|A| ^re fermionic complex fields. Therefore we can write 
where we have introduced the superfields <i>i,.-- ^|a| (* = 1, • ■ • , \M) 



(3.9) 



(3.10) 



These superfields can be seen as components of a supervector $ 

$1 



$ = 



|A| 



(3.11) 



|A| 



We adopted the conventions in the review by Mirlin We summarize super- 
symmetric formalism and notation in App. A. 
Now we can perform the average over H: 



(exp [-«<f>+i/<f>]) = exp 



(3.12) 



To convert this quartic interaction into a quadratic one we perform a Hubbard- 
Stratonovich transformation: 

J2 J., (<&+$, )(<i>+<f.) = Y.^AJ^JA, - B,J,,B, - 2P:j,,P,] (3.13) 



where 



■^i — Si, 



= XzXi, Pi = S*Xz, P^ = SiX* 



(3.14) 



and there is no sum over i. Note that Ai and Bi are commuting variables while 
Pi and P* are anticommuting ones. Now 



exp 
exp 



1 

+-B-JB 



(27r)- 



(27r)-V 



Vdet J 



VdetJ 



(3.15) 



p+ J^'- p-ip^ P~iP+ p 
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where a^, bi are real bosonic fields and pi is a complex fermionic field for any 
i = 1, |A|. Therefore 



(exp [--i$+i/<I>] ) = Jdadbdp*dpe i 



1 {a^ J^^a+b^ J^^b+2p+ J~'-p)-i'i>+ m 



where 



|A| 



Ri :— 



a^ p, 
Pt 



(3.16) 



(3.17) 



Ri is actually a supermatrix, containing both bosonic and fermionic variables. 
For such a matrix we can define the notion of transpose, complex conjugate, 
determinant and trace and it can be shown that the usual properties of the 
vector and matrix algebra hold. We summarize the notations in App. A. Using 
this formalism we have 



1 



n 



1 



E, ~ Rk J ^^^^^ SAet{E, - R,) 

(3.18) 



where 



Sdet(£;e " R^) 



1 



{Ee - a^) 

(E, - lb,) 

1 



1 * 
1 - Pi 



1 



^E^~RkJ^^ {E^ — Qk) 
Therefore 



1 - PtPk 



{E, - ai){E, ~ ib,) 
1 



{Es - ak){Es - ibk) 



(3.19) 
(3.20) 



(Gqq) ~ / da db dp* dp 



{Ee - ih) 



1 



{E, - ao) 



Pa Pa 



(E,-ao)(E,-iba) 



n[ 



1 - 



{E,-ai)(E,-ibi) 



(3.21) 



The integration over the fermionic fields can be performed exactly. Using 
the property: p| ^ (p*)^ = Vz we observe that 



1 - 



Pi Pi 



(E^-ai)(E,-ibi) 



exp 



Pi Pi 



{E,-ai){E,-ib,) 



(3.22) 



therefore the integration over p and p* reduces to the following expression 

-p'^[J-^+F{a,b)+F'{aa,bo}]p 



J dp* dp e 



det [J-^ ~ F{a,b) ~ F'{ao, bo)] (3.23) 



where F{ a, b) and i^'(ao,6o) are defined in (3^3.3). We obtain then the ex- 
pression (3.1). 
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4 Saddle point analysis 



In this section we shall deform the integral (3.1) over aj and bj so that they pass 
through certain complex saddle points. If we ignore the determinant in (4.9) 
and the kinetic term, we show that the resulting integrand has a double well 
structure, with the two wells of the same height. In the Sec. 5.1.2 (Lemma 6) 
we will see that the determinant actually suppresses one of the two saddles by 
a factor . 



Saddle points Observing the integrand in (3.1) we remark that the factor 
—W^A in a^J^^a + b'^J^^b forces the fields a and b to be approximately con- 
stant. Therefore if we ignore the determinant, the leading contribution to the 
integrand (hence also to the saddle structure) is then 



UE,-ib)] 


|A| 


[ [Ee - a) \ 




constant (a^ = a, 







|A| 



(4.1) 

b for all i) and we defined 

I 

h{b)="--\n{E-ib). (4.2) 

Note that in this approximation the saddle points for GUE and Random Band 
Matrix are the same. The critical points of /i and /2 are given by 



b. 



^ J- _l_ Z ^2 

~%^iy _l_ 



(4.3) 



where we defined 



Note that E satisfies 



E := Er - iti ■■= — - — 



E-£ = £\ ££* = l V|£;|<2. 



(4.4) 



(4.5) 



Spectrum Note that, if < 2 the saddle Os, 6s have non zero imaginary 
parts even as e | 0. For > 2 + 0{W~^) we expect that the density of states 
is smaller than any power of W^^, for W large. 



Contour deformation We deform the integration contour in order to pass 
through a saddle point. To avoid crossing the pole = Eg, we have to pass 
through the saddle Qs ~ £■ On the other hand the choice for bg is arbitrary, 



9 



as there is no pole in b, but it turns out (see Sec. 5.1) that bs = —i£ is the 
dominant contribution. Note that 



Hence the Hessian at this saddle point is 



Lemma 2 We perform inside (3.1) the translation 

aw ai + £ 



(4.6) 
(4.7) 

(4.8) 



bj bj — iS Vj' G A 



and take the limit e | 0. The integral can then be written as 

{G^o),=o = J dMa,b) det[l + (I? + I?[,)B] e^o'+^.^A (4.9) 

where the measure (i/iB(a,5) ^.l(\4y\ ) has covariance B given by l4-''')- The 
factor exp[^^ Vj] ((-11) is what remains in the exponential after the Hessian 
has been extracted, det[l + DB] corresponds to det[J~"'^ + F] after the 

normalization factor det B^^ has been extracted. Finally ex^pV^ and D'^ (i-H- 



4.14) o.f^ ih-e contributions from the observable. 

More precisely we define the measure as dfiB{o-,b) = d^Bio) d^isib) with 

dnB{a) = VdetB e-^(a^s-ia) ^ 



dfiB{b) = VdetTe 



(4.10) 



and B is the Hessian around the saddle, defined in (4-'') ■ The normalization 
factor for the measure has been extracted from the determinant. The interactions 
are given by Vj = Vj{aj) + Vj{bj) and Dij = SijDi, where 



V,ia,) = I dtil-t) 
V,{b,) = I dt{l-t) 



{£* 



ta,f 
2 i^b.f 



D, 





- F{a + £,b-i£)^^] ^ £^ 

r-l 

dt 



{£* - tib.Y 

1 

{£* - ai){£* - ibi) 
ibi 



(4.11) 
(4.12) 



_{£* - taiY{£* - iht) {£* - tai){£* - ib.tf 
Finally the contributions from the observable are given by Vq and D'q where 

V^ = -H£* -ao) (4.13) 



{D'o)^j ^-F'{ao+£,bo-i£k 



1 



{£* - ao)i£* - ibo) 



(4.14) 
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The proof is a straightforward change of variables and a reorganization of 
the resulting expression. Note that for any < 2 there is no pole in a as the 
factor E — £ — ai is always at a distance at least £i from zero. 

For the special value i? = 0, a singularity in 6^ = 1 seems to appear from 
the factor 1 / i{\ — bi) in the argument of the determinant. This is not a real 
singularity as there is the same factor in the numerator outside the determinant. 
Nevertheless to avoid additional technical problems we avoid £' = in the 
following. This is the reason why we chose 77 > in X in Theorem 1. 



Properties of the Hessian. The Hessian B ^ (4/7), which is the covariance 



of the Gaussian measure after the translation, has now a complex mass term: 



(1 -£^)^2i^l-—j+iE^l- — ml + im\ . (4.15) 

Note that for \E\ < 2 the real part is positive and this ensures the conver- 
gence of the integral. In the following, as we will need to treat in a different 
way the real and imaginary part of B~^, we introduce the real covariance C 

therefore 

B-^ ^ + iml . (4.17) 

Note that C is positive as a quadratic form and pointwise. In momentum space 
C is written as 

C(k) =^ ^ 7 (4.18) 

W^2j:f^,{l-cosh) + imr/wf 

h^27Tj^, n,=0,...,|A|i/'^-l . 

When |A| "f 00, ki becomes a continuum variable ki G [0,2tt]. The spatial decay 
depends on the dimension. In the particular case of d = 3 

The covariance B has the same expression as C, but with an imaginary term in 
the mass. It is easy to prove that B decays in the same way as C. 



11 



f|a) 

1. 

/ 

y 


\ 








f,(b) 

y 


a 









Ei 28i b 



Figure 1: behavior of F\[a) and ^2(6) 



Properties of the interaction After the translation the functions /i, /2 
introduced in (4.2) become 



/i(«)= + 



/2(&) 



(4.21) 



Note that, after the translation there also constant factors arising from /i and 
ji which cancel. In the following we will insert absolute values in the integral, 
in order to obtain our estimates. We then have to study the behavior of 



i^i(a) 



-/2(a) 



(4.22) 



It is easy to prove that for \E\ < 1.8 Fi{a) has only one maximum, in a = 0, of 
height 1 (see Fig. [^). Note that when 1.8 < < 2 zero is no longer the maxi- 
mum of Fi{a) and this is why we restrict E to T given by (2.5). Nevertheless, 
there is still a single saddle point so we expect that by suitable deformation of the 
contour we should be able to extend our r esuh to the interval \E\ < 2-0{W-'^). 

On the other hand, for any value of \E\ < 2, F2{b) has two maxima, which 
do correspond to the two saddles, one in 6 = and one inb = 2£i. Both maxima 
have height 1 (Fig. |l|) . We will see in the next section that the seco nd m aximum 
is suppressed by a factor e"^ from the determinant appearing in (4.£). 



5 Finite volume estimate 



We prove now Theorem 1 in a fixed cube A of side W , with G A. We prove the 
boundness of pa{E) in Theorem 2, then in Theorem 3 we prove the bounds on 
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the derivatives and on \pt^{E) — psc\- these bounds follow by the same technique 
used for the bounds on pa{E), with some slight modifications. 

Theorem 2 For A as above, there exists a value Wq such that for all W > Wq 
and for all E € X, where X is defined in (2^), the averaged density of states 
Pa{E) is bounded uniformly in W and A 



\PAm < K 



(5.1) 



Theorem 3 For A as above, there exists a value Wq such that for all W > Wq 



and for all E G X, where X is defined in (2.5), we have 



\d'EPA{E)\ < C„ yn<no{W) 



\pAiE) - psc{E)\ < 



K 



uniformly in A and W . 



(5.2) 
(5.3) 



5.1 Proof of Theorem 2 



Inserting the absolute values in the expression (4.£) we have 



(G+o),=ol ^ / \dl^B{a,b)\ \dct[l + iD + D',)B]\ 



(5.4) 



The absolute values of d^B and det[l + {D + D'q)B] are bounded through 
Lemma 3 and 4 respectively. 

Lemma 3 The total variation of the complex measure is bounded by 

\d^lBia,b)\ < dncia,b) (5.5) 



Proof The measure dpsio-j b) can be written as 

detB-i 



\dp.Bia,b)\ = 



detC- 



dpc{a, b) 



(5.6) 



where the determinants are the normalization factors for the two measures and 
can be written as 



detS" 



det C- 



det 1 + irriiCl 



(5.7) 



Note that, for any normal matrix A, with Tr < oo, the following inequality 
is true 



|det(l + A)| < le^'-^l e^^^^^ 



(5.8) 
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In our case A = i6C, therefore Ti A is imaginary, and the norm of the first 
exponential is one. The second exponent gives 



Tr A+A = mjTT ^ mf 
The bound in (|5.5|) then follows. 



1 



— — < 



■y- 



(5.9) 



Lemma 4 The determinant of I + {D + Dq)B is bounded by 
\det[l + {D + D',)B]\ < Ki^ eT'-(^+^o)B 



(5.10) 



Proof The proof is obtained by applying (5.S) and repeating the same argu- 
ments as in the Lemma above. Note that we applied sup^f, 1-0(0, 6) + I?o(a, b)\ = 
K for some constant K independent from W . ■ 



Applying the Lemmas above we have 



,Tr {D+D'a)B 



(5.11) 



where we bounded | exp(F')| = l'^^* - aol"^ < and V in defined in ( [IT^ ) 



Partitioning the domain of integration In order to distinguish small field 
and large field regions we partition the integration domain by inserting 



as follows 



k=l 



(5.12) 



dficia, b) 



JT,(D+D[,)B 



k=l 



n 



d^c{a, b) 



,TTiD+D'„)B 



where xi^''] is the characteristic function of the set l'^ and 

I^ = \a,b : |aj|,|6,-6,v| < -L Vj,/e Aand|6o| < -^1 

I^ = L,b : |a^-|,|6,-6y| < ^ Vj,/e Aand|6o-2£,| < ^1 (5.14) 
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^\a,b : hj £ IRVj and 3 j e A s.t. \aj \ > \ (5.15) 

L Wa J 

/* = la, 6 : \aj\ < G A and 3 £ A s.t. |6j - bj> \ > 

I W^s J 

/'^ = (a,6 : \a,\, \bj - bj,\ < Vj,/ G Aand|6o|, |6o - 2£,\ > 

Small field region The first two intervals correspond to the small field region. 
Ti is the leading contribution and corresponds to the case when all a fields and 
all b fields are near zero. In this case the interacting terms of the measure do 
not destroy the log convexity of the Gaussian d/ic, therefore we can apply a 
Brascamp-Lieb inequality jl^ ||l^ which states 

Brascamp-Lieb Inequality: Let 

dfinix) =: dxi...rfa;Ar— 5— e-^^(") (5.16) 

where x = (xi,...,xjv ) £ R^, H{x) is a positive function symmetric under 
X —X, and the partition function is 

Z{H) =: j dxi...dxNe-^"^''^ (5.17) 

Then if H" > C^^ > the following inequalities hold 

dfiH (x) \x^\" < J dfic (x) Ix^l"" n>0 (5.18) 

dfi„{x) e^^^^) < / d^ic{x) e^-^'^) (5.19) 



where dfic{x) is the free measure with covariance C, f is any vector in IR^, 
and {f,x) = J2t fi^i- 

The second term corresponds to the case when all the a fields are near 
zero and all the b fields are near the second saddle 2£i (see Fig. 0). In this 
case we bound the interaction (trace and Vj factors) by sup norm. The large 
contributions are now suppressed by a small exp[— M^] factor, from the trace 
bound. 

Large field region The last three intervals correspond to the large field re- 
gion. In all theses cases we bound the interaction terms (the trace and Vj) by 
sup norm in terms of quadratic and linear expressions in a and b. The large 
contributions from this bound are then compensated by the small probability 
factor (as the large field region is very unlikely). Note that the b field bounds 
are more delicate because of the double well structure (see Fig. p. 

Below we analyze the integration restricted to each interval. 
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5.1.1 Small field region: leading contribution Ti 

We consider the leading contribution Ti. In the region all the a fields and 
the b fields are near 0. We apply 



(5.20) 



(5.21) 



ReV{aj) < K ja^f if \aj\ « 1 
ReV{hj) < K |6jf if \bj\ « 1 

and we bounded the Tr{D + Dq)B applying 

l^.l < Wj\ + \b,\ 

\D'o\ < 1 

Therefore we can write 

Ti < / dA.c(a,6)e^'S.(l'^^l'+">^l')e^S.(l'^^l+l''^l)^"x[/' 



Now we insert the cubic and linear contributions in the measure by this 
definition 



H{a) =: a' C'^a 



Z{H) =: J dae-^"'-''^x[I^] 



Z{H 



(5.22) 

(5.23) 
(5.24) 

(5.25) 

(5.26) 



Therefore we actually have to estimate the normalization factor of the interact- 
ing measure This is done through Lemma 5 below. 



The same definitions hold for the b fields. Now we have 
'dA*c(a,6)e^S^(l'^^l'+l''^l')e^^.(l"^l+l''^l)'^"\[/i] 
where the free partition function is 



Lemma 5 With Z{H) and Zq defined by (5.2S) and (5.2t ) we have 

Z{H) < e°(w) Zo (5.27) 
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Proof Let H{t) defined as 



H{t){a) =: a^C~^a 



(5.28) 



interpolate between H and C . Note that on Ii 



In 



\Z{H{1))] 




[z{Hm\ 


-i 



< 



< K 



(5.29) 



where we defined dfijj(^fj{a) as in (5.24). In the last line we used Brascamp-Lieb 
( p.l8D together with 

H" > C-^ - fnij = C^i > (5.30) 
In general we will use this definition 



which is valid on Ii and ior f = O 

of C/ for / a constant < / < 1 or, when A is a set of cubes, / a diagonal 
matrix constant on each cube. Now 



z{Hm 

This ends the proof. 

Applying Lemma 5 we have 



< Zn 



(5.31) 



Ti < e°(i^) < K 



(5.32) 



5.1.2 Small field region: contribution from the second saddle 

In this section we show that T2 < e~'^^. This means that the fields have actually 
the same behavior as in a large field region. Note that this property holds only 
in three dimensions. 

In the interval all the a fields are near and the b fields are near the 
second saddle 2£j. RecaU that £i = y/l- E'^/A. Note that for aU \aj\ < W~i 
and \bj — 2£i\ < we have 



ReV{a,) <^faa 



Re V{bj) < ^ fbb^^ + (1 - fb)ml2£, {b, - £, 



2 •"--J 
2 



(5.33) 
(5.34) 
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with fa = fb = 0{W~^). Note that for b there is a hnear contribution coming 
from the translation to the second saddle. Moreover can be bounded Ti{D + 
Dq)B applying Lemma 6 below. 

Lemma 6 If \aj\ < and \bj — 2£i\ < W^i , then the real part of [{D + 

DQ)B]jj is bounded by 

Re {D + D'o)jBjj < -c W'^ (5.35) 
where c > is some constant independent from W . 



Proof Note that 

Re {DB).ij = [Re Dj Re Bjj ~ Im Dj Im B^j] (5.36) 
The key point is that, for aj near zero and bj near to the second saddle we have 

1 



Re ^ - ml + O 

1 



Im Dj = _ + O ( -i^ ) (5.37) 



Note that this estimates are not true in other regions. If both a.j and bj are 
near zero Dj ~ while for aj or bj far from the saddle we can only say that 
\Dj\ < const. For D'q we only need to know that for any oq, 6o G IR 

D', = 0(1) (5.38) 

Now, by simple Fourier space analysis we see that 

ReB,,>c^ (5.39) 



Inserting these estimates in (|5.36| ) the proof follows. ■ 
Inserting all this results in T2 we have 

T, < e-^W^ [ d^c:(a,5)e'#[^"S,a|+/.i:,&?]e(i-/.)™^2£,i:,(b.-£0 (5.41) 



We insert the quadratic terms in the measure. The normalization ratio are 
bounded using Lemma |^ below 

Lemma 7 For any < / < 1 we have 

det < ^ eOi^) (5.42) 



where Cf is defined in (5.3C). 
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Proof Diagonalizing the matrices we can write 

Y2j:U{l-cosk,)W^+m^^{l-f) 



1 o(^) 



< e 

- 1-/ 



where we defined C° as the covariance C/ where the zero mode has been ex- 
tracted. This ends the proof. ■ 



Therefore we can write 



<gO(w-^) e--W'< g-cw (5 44) 

where we inserted |A| ~ and fa = fb = 0{W^^) and we apphed 

'dA*c,,(6)e(i-^^)"'2^*^^(''^-^') = 1 . (5.45) 



5.1.3 Large field region 

This is the region selected by the intervals (one a fields large) (one pair 
of b fields with \bj — bjr \ large) and (all b fields far from both saddles). We 
apply the following inequalities 

Re {D + D'^),Bu < sup\{D + D^),| \Bu\ < O (^^^ (5.46) 
Re Via,) < "^faa] (5.47) 

9 

ReF(6,) < -^f,b] + 0[l-h) (5.48) 

with 1/2 < fa < 1, fb ~ 1 — W^^. These estimates are true for any value of 
aj and bj G H. On the other hand, when we are in the interval I^, all b fields 
must be far from both saddles the interaction in exponentially small, therefore 
we gain an additional small factor: 

ReV{b,)\j. < ^ + Oil -h)-c (^^) (5.49) 

Note that the factor Oil — fb) comes from the contribution of the second saddle 
(see Fig. |l|). Therefore we can write 

r3 + r4 + T5 < . (5.50) 
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We insert the quadratic terms in the measure: 

(1 - /a)(l - Jb) 

dfiCf,^ (a) dfiCf^ (b) {x[l'] + X[I^] + e-^^''^"*) (5.51) 



where we defined C and C as in ( 5.3C ) and we apphed Lemma ^ 



To bound the contributions from / and / we apply the following Lemma. 

Lemma 8 The probability of having one \aj\ > W~i or one pair \bj — bj'\ > 
is exponentially small 

dficAa) X[I^]< W^e-'^'^'^ (5.52) 
dl^cAb) X[I^] < iy«e-^^'^"* (5.53) 
Proof We consider first the integral for a 

f dfiCfM x[l'] < E / dl^CfM x{Wi\>W-i) 
J . J 

- X / '^l^Cf, (a) r e 



where we applied {Cf^)jj = 0{1/W'^) and we set x = 0(W~'sW'^). The same 
proof holds for the b field. In this case the presence of a difference bj — bj' is 
crucial to ensure that the factor [(C/Jjj + {Cf^)j'j' — 2{Cf^)jj'] is of order 
and does not depend on the mass (which could be very tiny for C/J. The factor 
comes form the sum over j and j' . m 

Putting together all the factors we have 

1 0(^1^^) o{i^+a-f.m. 



^^^^^^^^ ^ 



where we have inserted |A| = 0{W^), fa = 3/4 and /b = 1 - {\/W^). Note 
that there is an additional factor from the zero mode (1 — fb)^^ of the 
determinant. 
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5.1.4 Sum over the different regions 

Summing the bounds on different intervals we have finally 



< 



l + e-^^" +e 



< K (5.56) 



This completes the proof of (^.l|) . 
5.1.5 Large volume 

It is straightforward to extend the above estimates to the case when A is a union 
of cubes. 

Corollary 1 The density of states in a union of cubes A is bounded by 

Ipa(£^)I < e°(^) (5.57) 

Proof In each cube we apply the bounds above. The result is written as a 
quadratic form exp[?;"^C/w] where w is a vector which depends on the bounds on 
each particular cube and / is now a diagonal matrix which is constant on each 
cube. The key point is that 

where Aat is a Laplacian with Neumann boundary conditions on the cubes, and 
decouples the cubes automatically. Now we can perform the estimates in each 
cube separately. This completes the proof of (5.57). ■ 

5.2 Proof of Theorem 3 

To prove this result we integrate by parts to generate perturbative terms. To 
control the remainder we apply the bounds of Theorem 2. 

5.2.1 Semicircle law 

We prove that p{E) = psc with a precision of order W~^: 

1 



PAiE)^psc + 0{^] . (5.59) 



Note that pa{E) = — ilm(G[j'Q) therefore we have to study 

(G+ ) = / dpBia, b) — i ^e^ det[l + [D + D',)B] (5.60) 

J (c ~ ao) 
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We have to perform a few steps of perturbative expansion on the observable 
{£* — ao)~^ and Dq. These are more clear if we write the determinant as a 
fermionic integral. 

det[l + {D + D'a)B] = J dtiB{p\p) e-P'P'' (1 - D[,pyo) (5.61) 

where we defined 

dfiB{p*,p) = detS e""*-^"'" (5.62) 



*pD = 5]p>,D, 



(5.63) 



and Dj and Dq are introduced in (4.12) and (4.14). The density of states is 
then written as 



(G+) - / d^lBia,b,p\p) e^-"*"^ Oo 



(5.64) 



where we defined dpB{a,b, p* , p) = dpB{a,b)dpB{p* , p) and the observable Oq 
is 

1 



Oo 



{£* - ao) 



(1 - D'qP*po) 



(5.65) 



£ + ao dt 



1 



1 



{£* - taof {£* - ao) 



DqPoPo 



The first term is a constant and gives the semicircle law —-lm£ — psc- Note 
that we apply 



dpB{a,b,p*,p) e''-'''f"' = 1 
The remaining two terms give the corrections 



(5.66) 



6pi = / dpB{a,b,p* ,p) e 



V-p'pD 



ao / dt- 



lo " {£* - taof 
5p2 / dpBia,b,p*,p) e^"'' [-D'„pyo] 



(5.67) 
(5.68) 

Estimate of dp2 We first consider the estimate on the second integral, as 



it is the easiest one. We partition the integral domain inserting (5.12) as in 



Sec. 5.1 and we perform the fermionic integral in a different way depending on 
the region. 

Near the first saddle (interval I^) we apply 



dpBip*,p) e-"*"^ p*oPo = 



1 



B- 



D 



dct[l + I?B] 



(5.69) 



00 



22 



It is easy to see that for aj and bj near zero Dj ~ and \{B ^ + D)qq\ = 
0{W~'^). Therefore we have 



diiB{a,b) e^dct[l + DB] 



-D' 



i£*-ao) 



(5.70) 



where we apphed the same bounds as in Sec. 5.1.1. 

In the other regions /'^, fc 7^ we cannot apply ( 5.69| ) as {B^^ + D)^^ is not 
well defined (D is big and may cancel B^^). Therefore we apply 



d^Bip*,p) e'P'P^ p*aPo = detM 



where M is the matrix 1 + DB with the row substituted with B: 



M„- = (1 + DB)^J 



Mo, = 



DiBio 
Boj 



i ^ 



(5.71) 



(5.72) 



If we apply (5.8) we obtain the same bounds as in Sec. 5.1. Therefore performing 
the same bounds as in Sec. 5.1.2-5.1.3, we have 



J2\^P2il'')\< 



k^l 



Hence 



(5.73) 



(5.74) 



Estimate of S pi Now we consider the first error term. Before inserting the 
partition (5.12) and integrating over the fermionic integrals we have to perform 
one step of integration by parts 



5pi=^Bofc J dpLB{a,b,p*,p) 



(5.75) 



dt- 



1 



oak oao Jq (t* ~ taoy 



Note that \^Vk\<\ flfcP + lofcp. In the region around the first saddle (/^), 
applying the Brascamp-Lieb inequality ( 5.18 ), these fields give a factor W'^. 
In the other regions they are bounded by the exponential mass decay. The 
contribution from p^Pk is estimated as in Sp2 above. Therefore 



\Spi 



< 



J2\Bok\ 



\B 



00 



< 



(5.76) 



where we apphed J2k \^ok\ < const and |-Boo| = 0{W This ends the proof 
of Kb^ 
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5.2.2 Smoothness 



Now we consider the derivatives. Note that it is easier to compute the derivatives 
on the starting expression (GJq) than directly on the functional integral ( 5.60| ). 
The derivative at order n is given by 

91^PA(£;) = -n!(-l)"^Im((G")oo)«Im ^ (Goj, ...G,„o) (5.77) 

Applying the supersymmetric approach and the saddle point analysis as in 
Sec. 3-4, we can write for instance R{x) as 

R{x) = {GoxG^o) = (Co Ox) SUSY - C^oisusY {0'^)susY (5.78) 
where we defined 

{F(aAp\p))suSY - I diiB{a,b,p\p) e^'-p' F{a,b, p* , p) (5.79) 



and the observables are ( 5.65| ) and 



Ox = 



1 



{£* - ax) 



(1 - D',p*,px) 



(5.80) 



A similar formula holds for the general case. 

We perform now integration by parts starting from Oq until we have a path 
of connected vertices that connects to j or we have enough vertices to extract 
a factor for each observable Oj . This factor ensures that we can sum over 
the position of j inside the cube A. 

Note that, as in general we will have to estimate products of fields, both 
fermionic and bosonic, we will need the two Lemmas below. 

Lemma 9 Let consider the average of the product of p fermionic fields 

p 



diiB{p\p) e-"*"^ 



n 

.fc=i 



PikPi 



(5.81) 



Note that ik and jk are not necessarily equal. This integral gives different esti- 
mates depending on the region we are considering. If a and b are near zero we 
have 



dpB{p\p) e-" P^^xV' 



.k=l 



< 



pi 



\det{l + DB)\ x[l'] (5.82) 



On the other hand, in the other regions I" , s ^ \ we have 

p 



dpBip*,p) e 



n 

.k=l 



= a [detj/M]x[/^ 



(5.83) 
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where M is the matrix 1 + DB where the rows ii,...ip are substituted by the 
corresponding rows of B, and the columns ji,...jp are substituted by the corre- 
sponding columns of DB: 



Mij = (l + DB) 



D,B 



i ^ ii, ...ip, j ^ ji, ...jp 
i ^ ii, ...ip 



(5.84) 



Finally I are the set of indices I — {ii, ip} J — {ji, Jp} and detj/ M is 
the determinant of the matrix M without the rows Ji,,...jp and the columns 
ii,...,ip, and a is a sign. This new determinant can be bounded as usual. 



(5.85) 



Proof To obtain the first bound ( 5.82 ) we apply ( ^.69 ) for a product of 
fermionic fields. The result is the determinant of a p x p matrix whose ele- 
ments are {B~^ + -D)~-^ with i — ii, ...,ip and j — ji, ...jp. This determinant 
easily bounded by p\^sup\{B-^ + D):r.^\. Applying K^-^ + D)"-^] = 0{W- 
we obtain the result. 



IS 



The second expression (5.83) is easily obtained using the anticommuting 
properties of the fermionic fields. Finally (5.85) holds because the only error 
terms come from the absence of a term 1 in p diagonal elements. ■ 



Lemma 10 We consider the integral 

f ^ 
laini, ...,np) =: j d^xc{a,b) \e^\ |a 

k=l 

where p > 0, Uk > for all k and n — n^. Then 

7,(ni,...,np)[/i] < n\ (^^) 
/a(ni,...,np)[/2] < i^" e-^^ 



q>2 



(5.86) 



(5.87) 
(5.88) 
(5.89) 



// instead of a fields we have b fields the result is the same, but in the large field 
region we pay a larger factor, because we have a very small mass remaining in 
the covariance 



cW7 



q>2 



(5.90) 



25 



Proof As in Sec. 5 we partition the integration domain 1 = J2q=i xl^"^]- 
When we are near the first saddle ( /^) we write 

> < i ^ |a,J" (5.91) 



Now we can apply Brascamp-Lieb inequality as stated in (5.18). 

In the region near the second saddle ( we can bound the field a by a 
constant. 

In the large field region we bound the fields a using a fraction of the expo- 
nential decay of the mass term 

|a,r < V^e5*"'"? (5.92) 



.VsJ 

where 5 > is a small constant (5 < 1 which must be smaller than the mass 
of Cf^. Note that, for the b fields in the region P or I"^ S must be of order 
6 — 0{W~^) as this is the mass of C/^. This completes the proof. ■ 



6 Infinite volume limit 



In this section we shall establish bounds on (G^q) and the exponential decay 
of R{x) uniformly as A "f This is done by a standard method (see |^ 

or ||l^, ch.III.l) in statistical mechanics called the cluster expansion. These 
expansions are possible when there is a single dominant saddle point (in our 
case a — b = 0) whose fluctuations are close to that of a massive Gaussian i.e. 
a Gaussian whose covariance B has exponential decay. We are going to use a 
standard expansion with a few modifications. By supersymmetry some terms of 
the expansion are one (see Lemma 11) thus simplifying the expression. On the 
other hand for technical reasons the treatment of the covariance in the measure 
is slightly different form the usual one. 
We prove the following theorem 

Theorem 4 There exists Wq such that for all W > Wq hmA|^3 Pa{E) is 
bounded in X uniformly in W 



lim pk{E) 



< K 



for some constant independent from W . Moreover 

lim dlpK{E) < Cn yn<no{W) 
lim pa{E) = psc{E) + O {W-^) 
uniformly in A and W . In particular, for x ^ Q 



lim R{x)\ 

AT^3 



< e 

- 



(6.1) 

(6.2) 
(6.3) 

(6.4) 
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Outline of the proof: Note that the exponential decay (4.20) of B means 
that regions at a distance higher than W are approximately decoupled. As the 
observable depends only on the fields at zero or at zero and x we expect that 
all interactions take place in a volume of order around i = 0. To exploit 
this fact we partition A in cubes of side W forming the lattice C. For that we 
introduce the function 

Xa(«) = < ^ ^, . (6.5) 
I otherwise 

satisfying X]Ae£XA(*) = Xa(*)- the following we call root cube the cube 
containing i = and we denote it by Aq. 

The cluster expansion expresses (Gjo) and R{x) as a sum of finite volume 
contributions. Let 1" be a union of cubes A containing zero (zero and x when 
we estimate R{x)). Then the cluster expansion gives 

{Goo),=o - E ^y^y (6-6) 

where Zy is a functional integral over fields Oj, bj, pj and Pj, j G Y, and 
cy is a coefficient which is exponentially small in the length of the shortest 
tree spanning the cubes of Y. The main purpose of this section is to give a 



precise description of (6.6) and provide estimates to establish the convergence 
of the sum. Note that in conventional statistical mechanics there is usually an 
additional factor Zy with Y'^ = A\Y. However we show in Lemma 11 below 
that Zy = 1. 

The factors Zy are similar to a partition function except that there are 
derived vertices a^, &^ or p*pj{aj + bj) present. We shall show, using the ideas 
of Sec.5.1.5, that Z!y < (iCW-s)!^!. 

Lemma 11 If we restrict to the set of cubes Y'^ ~ h\Y and there is no observ- 
able contribution we have 

Zy. = J dpB^.ia,b,p*,p) eS.eA'(vS-P>.o.) = i (g.T) 



where dpByc {ci,b, p* , p) is defined after (5.64) o,nd By is the covariance B 
restricted to the volume Y'^. 



Proof We perform the translation aj aj—E, bj bj+i£ for all j E A'. Note 
that, for a general A' C A, B^/ ^ -W^A + (1 - £^). Therefore the translation 
gives some linear and constant terms. The constant terms are cancelled when 
we add the contributions from the a and b fields. By performing the inverse 
Hubbard-Stratonovich transformation we obtain 

= Q e**^(^-+^-^)*^ = 1 (6.8) 
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where the average ( )i is computed with the probabihty distribution (2.1) with 
covariance J instead of J, with Jy = B^^ + £^ . The matrix A is a diagonal 
matrix 



i k 

This completes the proof. 



(6.9) 



6.1 Cluster expansion 

We derive the cluster expansion formula. The result is stated in Lemma 12 
below. We construct the expansion by an inductive argument. 

First we want to test if there is any connection between the root cube Ag and 
some other A g A. For that purpose we introduce an interpolated covariance 
B{si) with < si < 1, which satisfies B{1) — B while B{0) decouples the root 
cube Aq from the rest of the volume. The easiest choice for B{s) is B{s)ij = sBij 
for i G Aq and j ^ A\Ao, or vice versa, and B{s)ij = Bij otherwise. For 
technical reasons we choose the following (less natural) interpolation rule 



Bisi)-' ^C{si)-' +imi (6.10) 



where 



> I siCij if i e Ao, j e A 7^ Ao, or vice versa -, -.x 

C{sihj = < „ . (6.11) 

I Cij otherwise . 

The reason we use this definition of B(s) is that we do not want to mix the real 
and imagin ary p art in B~^ in order to apply later the same estimates of Sec. 5. 
[) is equivalent to the definition 



Note that (|6T1 



C(si) = SiC + (1 - Si) [CaoAo + C^A^As] (6.12) 



{Caa')ij = ^ [XA(«)CjjXA'(j) + XA'(i)CyXA(j)] (6.13) 

where Ag = C\Aq. Therefore C(s) is still a positive operator, as it is a convex 
combination of the positive operators C and Caa- This fact is essential to 
ensure the convergence of the integrals. With the interpolated covariance we 
define 

Fa[si] = J dfiB(s^){a,b,p*,p) eS.^-=(^-''>^^^) Oo • (6.14) 

Note that for si = 1 Fa[si]sj=i = {Goo)^-q- Now we apply a first order Taylor 
formula to i^A[si] 

Fa[siU^i = Fa[si]si=o + / dsi ds,FA[si] . (6.15) 

Jo 
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The first term Fa[si]s^=o — Faq corresponds to decoupling Aq from the rest of 
the volume. The derivative in the second term of ( 6.15| ) gives 



(6.16) 



Using integration by parts we have 

J ds,dfiB{si)ia,b,p*,p) = j dpB{si)ia,b,p*,p) ^^^^^(si). 



_6__S_ _S_J_ J S_ 

Stti 5a j Shi Sbj Sp* Spi 



(6.17) 
(6.18) 



The derivative S/Sai may apply to exp[T^i — p*piDi] or to the observable Oq 
(this only for i — 0): 



" y Vi y ■ Vi 
Sai 



oai 



Sat 



—Di{aibi) 



{-p:pi)e~^^*^p'^P^''' 



Oo^{-p*oPo)^D'o 
dao dao 



(6.19) 
(6.20) 
(6.21) 



The same definitions hold for 6/dbi. The fermionic derivative 6/6pi may apply 
to exp[~p*piDi] or to the observable Oq (this only for i = 0): 



5 Pi 

Similar formulas hold for S/5p*. Therefore 



(6.22) 
(6.23) 

(6.24) 



where 

^A[si]((iiji)) = / dpB{si)ia,b,p*,p) 



6 6 



S S 



S 5 



Sa^ Saj^ Sbi^ Sbj^ 5p*^ Spi^ 



(6.25) 



Let us consider the propagator ds^B{si)i-^^j-^ extracted by the Taylor formula. If 
we choose the easiest interpolating rule, that is B{s)ij — sBij when i and j are 
in different cubes and i or j E Aq, the derivative is not zero only for ii E Aq and 
ji ^ Aq, or vice versa. Hence 9si5(si)iiji connects explicitly Aq to a different 
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cube. With the definition ( |6.10| ) the derivation is different and instead of one 
Hne we extract three 

ds,B{s,),,,, = ds, L.. . .^2 ) (6-26) 



C(si)-i+zm|, ^^^^ 



where we used ( |6.11 ) and G(s) is 

For each term (ki,k[) with fci G Ao and fcj^ S Ai we say there is a strong 
connection between Aq and Ai. We denote this by drawing a hne from Aq 
to Ai. Note that these points do not correspond to any derivative inside the 
functional integral, as the only derivatives occur on zi and ji . If ii and ji belong 
to some cube A ^ Aq U Ai they give some additional strong connections. 

Therefore the first step of the induction extracts a link li, associated to 
four points ii, ji, fci, fci, connecting Aq to a set of one, two or three new cubes 
depending on the positions of ii and ji. We call this set the generalized cube 
Ai. The different possible links inside Ai are shown in FigJ|. Now 

ds,FA[si] = Gisi),,k^Ck,k'Gis,)k'^^^ F^[sl]ii^^n)) ■ (6.28) 

ihdi) 

Note that the fimctional integral after ds^Bs^ has been extracted is function 
only of (ii, ji) and not of (/ci, k'l), as only the first two indices correspond to a 
functional derivative inside the integral. 

Now we fix the points Ji), (fci, k'l) corresponding to a strong connection 
between Aq and Ai. We want to test if there is any connection between the set 
Ao,i = Ao U Ai and any other cube A' G A\(Ao,i). We introduce then a new 
parameter S2 in C(si) in the functional integral: 

^, . ]S2C{si)ij if i e Ao.i, j ^ Ao.i, or vice versa 
G(si,S2)y = < „, , ^, . (6.29) 

I L[si)ij otherwise . 

As for C(si) we can write C{si, S2) as a convex combination of positive operators 

C{suS2) = S2C(si) + (1 - S2) [GA„^Aoa(^i) + ^Ki^M.S'^ (6-30) 

where C(si) and Caa(si) are positive. Therefore G(si,S2) is still positive. 
Then Ff^[si\{{iiji)) = Fa[si, S2]((iiii))s2=i- We apply again the first order Taylor 
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formula 



-FA[si,S2]((nil))s2 = l = Fa[si, S2]i{iljl))s2=0 + / ds2 9s2i^A[si,S2]((»lil)) 



(6.31) 

As before Fa[si, S2] ((nii))s2=o corresponds to a functional integral restricted to 
Ao,i: -Fa[si, S2]((iiji))s2=o = F^^^[si]{(hji)). The other term gives 

ds2FA[si,s2]{{iin)) = ^ [G{si,s2)t2k2C{si)k2k;^G{si,s2)k'^j2] 

(12 J2)(A:2,A:2) 
fc2eAo,i,fe2eAo 1 

X FA[Sl,S2]{{^ln)A^2j2)) ■ (6.32) 

We repeat this argument until we construct all the possible connected compo- 
nents containing the root cube. This is a finite sum, for A fixed. 



Definitions We give now some more precise statements. We define a gener- 
alized cube A as a set of one, two or three disjoint cubes in A. A generalized 
polymer Y is then a disjoint set of generalized cubes A. A tree T on F is a 
set of links ^i,.. connecting the generalized cubes in Y and forming no loops 
(See Fig.||). We call the set of all these cubes the polymer Y contained in Y. 
Each link Ir corresponds to four connected vertices v, jV, kr, k'^. The corre- 
sponding propagators are Gi^^k^, Gk^,k'^ and Gk'^.j^ defined in ( 6.34| ) and the 
corresponding links are shown in (see Fig.|^) . Note that the same A and tree T 
may correspond to several different polymers Y (see Figj^). 
We apply then the following formula. 



Lemma 12 {G'^q) ^_ „ can he written as 



(G^o),.o- E E n f<^'r E ^t{s) 
YBAo Tony r=l ii^-.l,-') 



(6.33) 



\Y\-1 



n G{s),^k^Ck^k'G{s)k 



r=l 



FY[s]{{ir,jr}) 



where we sum over the generalized polymers Y and over the ordered trees T on 
Y with root Aq. Note that the ordering on a tree is the ordering on its links 
Ir, r = l,...|y| — 1. The product over r is then the product over the links in 
the tree. The points {ir,jr) {kr,k'^) fix the number of cubes inside each A €Y 
and the links connecting them. Each tree link Ir is associated to the parameter 
< Sj. < 1- The product Gi^krGk^k' Gk' corresponds to ds^B{sr) for each 
link Ir (z T and the factor Mt{s) is the product of s factors extracted by the 
derivatives drB{sr). Finally JV[s]({ir, jr}) *s the functional integral remaining 
after the propagators ds^B{sr) have been extracted. 
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Figure 2: The two polymers and tree structures in b. and c. both correspond 
to the same generaUzed polymer Y and tree T in a. 



More precisely, for each r the points {ir,jr) {kr,k[) must satisfy the con- 
straint kr e Ar' for some r' < r, k'^ € A^, and ir,jr ^(^y belong to any A^.' 
with r' < r. The propagators C{s) G{s) are defined as 



G{s) 



" 1 + im^C{s) ' 



where 



1 
1 





if 3 r s.t. i,j G Ar 

if 3r' < r s.t. i e A^, j € A^ 

if i e F, j^Y . 



The remaining functional integral is 

FY[s]i{ir,jr}) = J dfXB{s)ia,b,p*,p) 
S 5 S 5 6 6 



n 



6ai^ 6aj^ 6bi^ 6bj^ 6p* 6pi. 



where B{s) is 



B{s)-' := C{s)-' 



(6.34) 



(6.35) 



(6.36) 



(6.37) 



Note that, as we constructed the tree, the order on the tree lines ensures that, 
for each 1 <r < n, the tree line Ir connects Aj. to some Aj.' with r' < r. 
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^0 "1 

a. 



Figure 3: some examples of links of type 1,2 and 3; the two point vertex is a 
filled dot, the three point one is an empty dot 



Proof The proof follows directly from the inductive procedure we explained 
above. ■ 



Now we can bound [(Goo)! inserting the absolute value inside ( |6.33|) . There- 
fore we have to bound 



Y TonY «=1 



{K,k'^). 



\Mt{s)\ 



(6.38) 



'\Y\-1 

n \g{sKk\\Cu^k\\g{^)k^A 



6.2 Bound on the functional integral 



\FY[s]{{^r,Jr})\ ■ 



We bound the remaining functional integral |^V[s]({ir, jr})| by a generalization 
of the arguments for the finite volume case (Sec. 5). 



Definitions We call Vd the set of vertices derived by the cluster expansion 
appearing in ( 6.3(\ ) : 



{j eY : 3 1<r <\Y\~l s.t. j = > or j = v} 



(6.39) 



For each j G Vd we call dj{a), dj{b), dj{p), dj{p*) the number of derivatives 
6/6aj, 6/6bj, S/5pj or d/Sp* respectively. As we see from ( |6.19 - |6.22 ) these 
derivatives apply to V{aj), V{bj) or Djp*pj. We also have a contribution from 
the observable Oq, but only for j = 0. 



33 



Note that, by the properties of anticommuting variables, dj{p), dj{p*) can 
be only or 1. On the other hand there is no limit to dj{a), dj{b). We call 
dj = dj{a) + dj{b) + dj{p) + dj{p*) the total number of derivatives in j. This 
number is actually fixed by the choice of {{ir,jr)}r- For each j G Vd we have 
to study 

Qd,ia)Qd,ib)gd,ip')gd,(p) (^gy,(a)+y,(6)(i _ p*p.Dj)^ (6.40) 

In the small field region we need to extract some structure. We can actually 
write the derivative as 

^ af'^ bf'^ pf"^ C{a^)C{bi) (6.41) 

'"j(p)>'"j(P*) 

where rj{a), rj{b), rj{p), rj{p*) are the number of fields remaining after the 

derivatives have been performed. Note that rj{p) and rj{p*) can take only the 
values and 1. On the other hand it is easy to see that rj{a) < 3dj{b), and the 
same holds for b. Moreover the parameter nj = dj + rj > 3, except for j = 
(and j = a; if we are considering R{x)). The factors C(a) and C{b) no longer 
depend on the fermionic fields. By analytic tools we can show that in the small 
field region 

\C{a)\ < K''^ djiay. |C(6)| < K'^^ dj{b)\ (6.42) 

In the large field region (that means in for q > 1), we do not need to extract 
the whole structure, as the fields are large and the small factors come from the 
probability. We only need to extract explicitly the fermionic fields therefore we 
write the derivatives as 

^ p;.(p) py.ip') C{aj)C{b,) (6.43) 

Again by analytic tools we can show that in the 

\C{a)C{b)\ < K'^' dj{ay.dj{by. e^^(")+^^('') (6.44) 

In this region we have no factors rj . Nevertheless wc define rj = if > 3 and 
r j = 3 — dj otherwise. In this way we ensure nj > 3 for all j G Vd- 
With these definitions and results we prove the following Theorem. 



Theorem 5 

\FY[s]{{ir,jr})\ < i^r' E n 



(6.45) 



where Ki, are constants and ta = X^jeVdnA'^j' "'A = X^jeV^nA' 
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Proof The proof is a generalization of the arguments for Theorem 2 and 3 in 
Sec. 5. Note that Lemma 3 and 4 hold also after substituting B and C with B{s) 
and C(s). Then we can introduce the partition ( 5.13| ) in each cube separately 



n 



(6.46) 



We perform the bounds in each cube as in Sec. 5. Note that, as the derivations 
bring several bosonic and fermionic fields out of the exponential we have to use 
some of the ideas of Theorem 3. 

First we have to perform the integral over the fermionic variables extracting 
the correct factors W^^. Now we ha ve ma ny different cubes, each in a different 
interval, therefore we cannot apply ( 5.82| ) as in Lemma 9, as {B^^ + D)^^ is 
not be well defined. We apply then a generalization of (( ^.83 ) which allows to 
extract also the small factors. 

We have to compute the integral 



dtJ.B{s}{p* , P) 



(6.47) 



n 



ieVd :rj(p) = l 



n 



jeVd :rj(p*) = l 
. rj(p)=0 



n 



pjpj 

rj{p}=rj{p*) = l 



We partition now the set of j G F as 



U - {jeY\j^v4 

V = {jeY\jeVd, r,{p*)^ lor r,{p) = l} 



(6.48) 



Note that the points j G Vd with rj{p*) = r.j{p) = 1 do not give any contribu- 
tion to the integral. For each cube we introduce r^ip) = X^iev^nA ^^'^ 
TAip*) = SiGVdHA Note that the fields pj are columns in the resulting 

matrix and the fields p* are rows With these definitions we have the following 
Lemma. 



Lemma 13 By multi-linearity of the determinant, the fermionic integral ( 6. 47 } 
above can he written as 



n 



rA(p)''^<PVA(p*)''^^P'^ 

|yA(p)+rA(p*) 



a det M 



(6.49) 



where a is a sign depending on how we order the rows and column in M , and 
M is the block matrix 



M 



Muu Mud 
Mdu Mdd 



(6.50) 
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Muu is 0, matrix corresponds to the elements {M)ij with both i,j G lA, M^d to 
the elements {M)ij with both i,i ^T) and Mud o,nd Mdu cire the mixed terms. 
Note that, if j e Vd and rj{p) — 1 but rj{p*) — j appears only in a column of 
Mdd and the element Mjj is not present in the matrix. Therefore we order the 
lines and columns of Mdd'- {Mdd)ij = {^)iiCj where U and Cj G Vd. With these 
definitions the matrix elements are 

w 



{Muu\^ = [l + {D + D'^)B{s)]^^ {Mudh^ - D,B{s),,^^^^ 

(6.51) 



where A(i) is the cube containing the vertex i. 

Proof The proof follows from the properties of the anticommuting variables 
and determinants. ■ 

Now we can insert absolute values inside the bosonic integral. If we bound 



det M as in (5^) we obtain the same bound as in Lemma 4, with an additional 



error term. The precise statement is given in Lemma (14) below 



Lemma 14 The determinant of the matrix M defined as in (6.5C-6.51) satisfies 
the bound 

IdetMl < ei^^^(^""-i) eOC^I+l^''') (6.52) 



Proof Using ( |5.8| ) we have 

IdetAfI < eR-Tr(M-l) gKTr(M-l)*(M-l) 



(6.53) 



Now applying the definitions (6.5C) we can write 



Re Tr (M - 1) = Re Tr (M„„ - 1) + Re Tr {Mdd - 1) 

Tr (M - 1)* (M - 1) = Tr (M^ - l)*(Af„„ - 1) + Tr [Mdd - l)*(Af<jrf - 1) + 
Tr [{Mdu - l)*(M„d - 1) + {Mud ~ \r{Mdu - 1)] (6.54) 

By inserting the definitions ( |6.51 ) and using the decay of Sy , it is easy to 
see that all terms are bounded by a constant per cube except for 



Tr(Afdrf - \r{Mdd - 1) < 51 ''a(p) + r^{p*) 



(6.55) 



This completes the proof. ■ 

Now we perform the estimates as in Sec. 5. As in Sec. 5. 1.5, in order to 
decouple the estimates on different cubes we write all the bounds in terms of 
quadratic forms like v'^Cf{s)v where v is some vector and Cf{s) = C{s)~^ — 
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/m^ > as in (5.30) but / is now a diagonal matrix which is constant on each 
cube. Then we can apply 

C(^) < , (6-56) 



N 



where Ajv is the discrete Laplacian with Neumann boundary conditions on the 
cubes. This operator decouples automatically different cubes. Note that to 



estimate the contribution from the small field region we have to apply (5.87) 
Lemma 10. To extract the factors W~'^^ in the large field cubes we use the 
exponential factors e~^. 



6.3 Sum over the clusters 

We perform now the sum over the clusters Y. We split the sum in several pieces. 
First fixing the cubes A we sum over the points jr, kr, k'^ in the cubes. Note 
that after this operation is done, there is still a small factor associated to each 
cube. The factorials arising from the combinatorics are beaten by a piece from 
the decay of GCG. The remaining piece of the decay is used to sum over the 
cube positions, following the tree structure. Finally we sum over the tree choice 
T using the fact that we have a small factor per cube. 

To perform all these bounds we need now to study the spatial decay of the 
propagators C, B and G. We know already the spatial decay of C(s) (see ( 4.191 ). 
The decay of G is given by the following Lemma. 

Lemma 15 The propagator G = (1 + imfG {s))^^ decays as 



|Gy I < Sij 



2 |nii 



-/- 



il 



where f — inf[l/2, and g is some constant independent from W . 

Proof By a Combes-Thomas argument we prove that 

1 



(6.57) 



R 



-1 



< 2 



(6.58) 



l + imfC{s) 

for R a multiplication operator defined as R\x exp[/ix] \x > and jl any vector 
with \ < gnir/W. Now G can be written as 

1 



dij — im 



'■C{s),, - mt (c{s) . ^ C{s)\ (6.59) 

V i + imic{s) 



,l + zmfG(s),,^ V ^ . /.J 

We need to study the decay of the last term. Actually we see that 



G{s 



:C{s) 



/i|x-j7| 



1 + imfC{s) 
(R-^C{s)R) . , 



R-^C{s) 



1 + imfC{s) 
r) {R-^C{s)R) 



C{s)R 



\iV,AW)\<\\V\\\\W\\\\A\\ < ^ 



(6.60) 
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where wc defined Vk = {R-'^C{s)R)ik, Wj = {R-'^ C (s) R) kj and A = + 
imfC{s))R. We chose jl such that i^\x — y\ = fl{x — y), and = |/2| < gnir/W, 
and fi < mr/2W. The last condition must ensure that the exponential decay of 
C{s) controls the exponentials from R~^ and R. This completes the proof. ■ 

Note that this Lemma gives also the decay of B{s) 

\B{sU < l^TT^e--l-l + ^Me-/^K-.I (6.61) 
as B{s) = {G-l)i/ml 



6.3.1 Extracting small factors 

Before performing the estimates we extract some factors from the propagators 
GCG for each tree line, to offset the factorials eventually arising in the estimates 
and to ensure we have a small factor per vertex. 



Factorials. Constant powers of factorials such as d^, for p fixed, can be beaten 
using a piece of the decay of CGC. Note that each tree line Ir connects different 
cubes, therefore we have cIa disjoint cubes hooked to the cube A by the tree 
T. When cIa is large, since we are in a finite dimensional space, many of these 
cubes must be very far from A. It is easy to see that, half of the d-A cubes must 
be at a distance from A of order Wd^. Therefore we gain a factor exp[— cd^] 
which can beat any constant power of factorials. 



|y|-i 



< 



n 



UA'- 



IP 



(6.62) 



Note that the constant K depends on p. 



W factors We need a factor W ^ for each field hooked to a derived vertex 
j G Vd- We extract then a factor from each G propagator: 



y|-i 



n \Gis)irkA\G{s)k'^j 



n(^ 



n \WG{s)i^kJ\WGis)k'^ 



n 



(6.63) 



Note that rij > 3 for all j £ Vd except for j = 0. As the each cube has volume 
W^, this ensure that we can choose the position of each vertex without paying 
any factor W. 
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Finally, after extracting a fraction e of the exponential decay and the factors 
we separate in the remaining factors the polynomial and exponential decay 



\Y\-1 



'r-jr 



n 

.r=l 



(6.64) 



where 



Gij 



j_ 

w 



Cij 



(6.65) 



where /' = /m^ — e is the remaining mass and d{A, A') is the distance between 
the center of the cube A and A'. 



6.3.2 Sum over the vertex positions 

Now using the decay — j\ in G and C we sum over the positions of all 
vertices inside their cube (the cube is fixed). Each line of the cluster expansion 

corresponds to four vertices ir,jr,kr,k'j., where kr and k'^. correspond to two 
point vertices and must belong to different cubes while ir, jr may belong to the 
same cube. For each j = ir or j — jr we distinguish two cases 

• j contracts to j' and j' has never been extracted before in the cluster 
expansion. Then we say f is new with respect to j. 

• j contracts to j' and j' has already been extracted before by the cluster 
expansion. Then we say j' is old with respect to j. 

We consider first the case kr 7^ ir and k'r ^ jr so that the factors Si^k,. and 
6j^k'r disappear. Note that we sum over the position of ir {jr) only when it is 
new. We consider the different cases. 

ir ^ jr and both ir and jr new Then we sum over ir and jr- 

E E G^rk. E ^'^'■K E ^k'rJr = {w i) {w ^) 

(6.66) 

Therefore we pay a factor W^^"^ for v and the same factor for jr- 
ir jr and v new, jr old Then we sum only over ir. 

E ^Kj^ E ^''^'K E Gi^kr<W^ (6.67) 

fc;GA(j^^) fcrGA(fc^) i,-eA(i^) 

Therefore we pay a factor for ir and no factor for jr- 
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ir ^ jr and v and old Then v and jr are both fixed. 

G.^k.W-^ < 0{1) 



(6.68) 



4 



where we bounded Ck^^k' < W 
ir = jr and ir new Then we sum over v- 




E 



\ir-kr\<\iT-k'^\ 



(6.69) 



Therefore we pay a factor for v. 
= jr and v old Then v is fixed. 

51 G.^fc.W^"' E ^Ki. ^ 0(1) (6.70) 

where we bounded Ck^^y < W~^. When v — or jV = it is easy to see 
that the same estimations hold. 

Note that for each j = v (j = jr), with j 7^ 0, we pay a factor VF^/^ when 
it is new and some constant K in any other case. Therefore, applying Uj > 3, 
we have 



K 



< 



K 



(6.71) 



Therefore we have a factor W~'s for each tree line hooked to j. This means we 
have a factor W~'s for each generalized cube A is Y . The case j = is special 
as rij > 0. Nevertheless the position of is fixed so that we do not pay the 
factor W^^^. Therefore for j = we have 



1 



(6.72) 



6.3.3 Combinatorial bounds 

The combinatoric inside each cube costs a factor 



(6.73) 



The factorials are beaten by a piece of the exponential decay of the tree lines 
( |6.62 ), while the constant factor will be bounded later when we will sum over 
the tree structure. 
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6.3.4 Sum over the cube positions and the tree structure 

We sum over the cube positions using the exponential decay and following the 
tree from the leaves towards the root. The result is K^"^^. The remaining sum 
is now 



\Y\ T orders K^Y 

unordered 



K 



(6.74) 



where we have split the sum over ordered rooted trees as the sum over unordered 
rooted trees, and the sum over orders. This last sum is performed by the integral 
over the interpolating factors (see |16 or jl7j, Lemma III. 1.1) We give here a 
sketch of the proof. 

Lemma 16 The sum over all the orders on the tree T is bounded using the 
interpolating factors Si as follows 



J2 U J ds,\MTis)\^l 

orders i—1 ^ 



(6.75) 



Proof We introduce the variables for all Aj G Y. Then we introduce 
the function F{e) =: n(ij)GT(l + where T is unordered. Now we perform 
the tree expansion as we did in Lemma 12 Sec. 6.1. We define £y (si) = siSij if i 
or J = and £jj(si) = £y otherwise. We apply the first order Taylor expansion 
and we go on until we extract all the Sij. The term proportional to all Eij is 
then 



n E n / 

{ij)£T orders i=l " 



(6.76) 



If we expand F(e) in powers of e we see that the term Jl(y")eT coefficient 
1. Therefore by comparing powers of e we obtain ( |6.75| ). ■ 

Finally the sum over the structure can be written as 



E E 

15-1 



T 

unordered 



E U 

dr >0 i=l 



E n^E 



(6.77) 



where we defined g — K/Wb < 1 and for each generalized cube we sum over 
the coordination number. We start summing from the leaves going towards the 
root. The leaves give 



< 02 



5^ 



if 52 



d>0 



1 

~9 



< 1 



(6.78) 
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The following step gives 



'I]5'=5-^<5^ if5^-^<l (6.79) 



Then we repeat inductively. Finally 



E E 9^^^<K + l (6.80) 



lYI T 

unordered 



where the constant conies from the bound for n = 0. This completes the proof 
of the first part of Theorem 4, namely the boundness of the density of states 

6.4 Smoothness and exponential decay 



To bound the deri vatives of the density of states (|6.2| ), and in particular the 
decay of R{x) (|6.4[ ), we perform the cluster expansion as in Sec. 6.1. 

For R{x), note that contributions where Y does not contain both and x 
are cancelled. When and x both belong to Y we can extract the exponential 
decay directly from the tree lines GCG. For derivatives d^pA{E) the idea is 
the same. Only contributions from Y containing all the observables are not 
cancelled. The fine structure (the factor W^^ in R{x)) are then extracted by 
a few steps of perturbative expansion as in Sec. 5.2.2. In the same way we can 
prove the semicircle law behavior (|6.3|). This completes the proof of Theorem 4. 



Appendix A 

Supersymmetric formalism 

We summarize the conventions and notations we adopted in this work (they 
are based on the review by Mirlin §). 

Grassmann variables A set of Grassmann variables and their complex con- 
jugates xi: Xi5 ■■■Xn Xn following properties: 

XtXj = -XjXi, X*Xj = -XJX^, X*X*j = -XjX* (A.i) 

(XD* = -X», ix^Xj)* = x:x*j (A.2) 

Jdx^^^J dx* 1 = 0, Jdx^X^^J dx* X* - ^ (A.3) 
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With these definitions we introduce a vector and its adjoint as usual 



X = 



= {X*i,--- ,Xn) 



(A.4) 



V XJV / 

Now x'^X is a real commuting variable and 
for any matrix M . 

Supervectors and supermatrices A supervector is defined as 
/ ^1 \ 



(A.5) 



^+ = {Sl,--- ,S*i„xl,--- ,Xn) 



(A.6) 



Xi 

\XN ) 

where Si are the commuting and Xi are the anticommuting components. Simi- 
larly a supermatrix is a matrix with both commuting and anticommuting entries 



M 



a (7 
P b 



(A.7) 



where a and b arc ordinary matrices while a and p have anticommuting ele- 
ments. We identify the element of a supermatrix by four indices M"j^^ where 
a,p specify in which sector we are: (0,0) corresponds to a (boson-boson); (1,1) 
corresponds to b (fermion-fermion) ; (0,1) corresponds to a (boson- fermion); 
(1,0) corresponds to p (fermion-boson). (i.j) identify the matrix element inside 
each sector. For example M^P^ = ay. 

The notions equivalent to trace and determinant are supertrace and superde- 
terminant 

StrM = Tro - Trfe, SdetM = det(a - o-fe" V) detb'^ (A.8) 

With these definitions we have 



/ 



StrlnM = lnSdetM 



d^*d^ e-*^^* = SdetM-i 



(A.9) 



(A.10) 
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J e-*+^^* = iM-^)l^ SdetM-i (A.ll) 



Note that some properties are different from that of the usual matrices, in 
particular: 



Sdet zM = SdetM 



(A.12) 



for any complex number z. 

Finally from these formulas one can derive the inverse of the supermatrix 
M ( [TtI ): 

^ [ -b-'p{a-ab-'py' b-'\l+p{a-ab-'p)-'ab-'^ ' ^ 
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